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Bosonic cascades formed by lattices of equidistant energy levels sustaining radiative transitions 
between nearest layers are promising for the generation of coherent terahertz radiation. We show 
how, also for the light emitted by the condensates in the visible range, they introduce new regimes 
of emission. Namely, the quantum statistics of bosonic cascades exhibit super-bunching plateaus. 

This demonstrates further potentialities of bosonic cascade lasers for the engineering of quantum 
properties of light useful for imaging applications. 

PACS numbers: 78.67.-n,71.35.-y,42.50.Ar,78.45.+h 

Introduction. —Terahertz frequency radiation is a valu¬ 
able resource for a number of imaging applications in 
medicine, security screening, and different scientific dis¬ 
ciplines [1]. For this reason a number of solid-state sys¬ 
tems have been considered as terahertz sources, where 
one typically aims to convert an optical frequency pho¬ 
ton into a terahertz photon. Even if materials with ap¬ 
propriate transitions can be found, such a procedure is 
inherently inefficient as the majority of the energy of each 
optical frequency photon is given up in favour of gener¬ 
ating a low-energy terahertz photon. Quantum cascade 
lasers circumvent this problem, where an optical quan¬ 
tum of energy can be used to generate multiple tera¬ 
hertz photons i®. While fermionic (electronic) cas¬ 
cades have been realised over 20 years ago Q, bosonic 
cascades were only recently designed theoretically based 
on excitonic transitions in parabolic quantum wells, pos¬ 
sibly being placed inside semiconductor microcavities Q. 

When using also a THz cavity, bosonic cascades offer the 
prospects of double stimulation of the emission: by the 
final exciton state population [b], Q and by the terahertz 
cavity mode population [H, @|. 

Stimulated transitions among macroscopically occu¬ 
pied bosonic quantum states are in the heart of bosonic 
cascade lasers (BCLs), which also act as polariton lasers 
spontaneously emitting coherent light in the optical fre¬ 
quency range 0 - Polariton lasers bring the advantage 
of ultra-low threshold power 00 and as versatile all- 
optical platforms 0 they are expected to constitute 
building blocks of future optical integrated circuits 0 , 
logic elements 0 or polarisation switches 0 - The 
quantum optical properties of polariton lasers have been 


studied theoretically 0 and experimentally jl9]. The 
quantum theory of polariton based terahertz sources has 
been addressed recently [20]. The quantum physics of 
bosonic cascades is yet to be explored. 

In this work, we open a new dimension for BCLs by 
showing that their emission produces a marked super- 
bunching, thanks to their specificity of making coexist 
multiple macroscopic coherent states. This makes for a 
difficult problem to describe quantum mechanically, that 
we can tackle here with stochastic quantum Langevin 
equations. This theory of bosonic cascades to describe 
their statistical properties not only confirms their depar¬ 
tures from other mechanisms of lasing 0, it also sug¬ 
gests that they could power imaging applications with 
resolutions much higher than ever conceived before [22| . 
thanks to combinations of several assets to go beyond the 
diffraction limit Him. 

Theory .—We consider a bosonic cascade with N en¬ 
ergy levels labelled from A = 1 to TV, as illustrated in 
Fig. HJ The parabolic quantum well makes the levels 
equidistant in energy [25j. We account for the dominat¬ 
ing transitions which are between neighbouring excitonic 
levels A +1 and A and assume for simplicity that all levels 
decay radiatively with the same rate, r _1 . We assume 
that the occupation of the terahertz mode fed by the in¬ 
terlevel transitions in the cascade is zero, i.e., there is no 
terahertz cavity in our model structure. This corresponds 
to the original version of BCLjEj], where the stimulation 
of terahertz radiation is achieved due to the macroscopic 
exciton populations in the bosonic cascade. This simpler 
configuration already manifests the non-trivial quantum 
effects of bosonic cascades. The system is thus described 
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FIG. 1. Illustration of a bosonic cascade, a) Structure 
schematic: a parabolic quantum well is excited with an op¬ 
tical pulse and emits light at a range of frequencies, b) En¬ 
ergy level diagram: the parabolic trapping potential engineers 
equidistant energy levels and transitions occur between neigh¬ 
bouring levels with transition element V. Bosons in each level 
can decay radiatively, with characteristic decay time r. 


by the density matrix p subject to a quantum Liouville 
equation: 


% = 27 £ ([«a^«a] + KMa]) 

A=1 

v N ~ 1 

+ j ([«Ua+i>M+a] + K«a + F’4+a]) . 

A=1 

(i) 

where a\ annihilates a boson from level A and V is 
the scattering rate [26]. Dotting this equation with 
Fock states and neglecting correlations between levels 
in the Born-Markov approximation leads to the Boltz¬ 
mann equations, which describe well the average popu¬ 
lations of each state and thereby provide the mean-field 
theory of bosonic cascades [b]. To go beyond this semi- 
classical description, and in particular to compute the 
photon statistics, correlations between states must be re¬ 
tained. A brute-force approach, however, is restricted to 
few levels and each with a modest occupancy. To access 
the most general cases of bosonic cascades, with macro¬ 
scopic occupation of the states, we therefore expand the 
density matrix on the natural basis for this problem, that 
of coherent states [27]: 


P= j 'P(a 1 ,p 1 ,...,a N ,p N )A(a 1 ,p 1 ,...,a N ,l3 N )dii, 

( 2 ) 

where: 


A(aq,/3i,..., olniPn) 


l a l, • • • ,«iv) (/?!*,. • - ,Pn\ /o') 


with V the positive-P distribution, which differs from 
the Glauber-Sudarshan distribution in allowing for non¬ 
diagonal coherent state projectors. The complex num¬ 
bers a\ and /3\ are independent variables covering the 


whole complex plane. The integration measure is dp = 
d 2 aid 2 /3i ... d 2 ajyd 2 /3N and for ease of notation we write 
a = (aq, /?i,..., ajv? Pn)- a ? n will be used to refer to an 
element of this vector (of length 27V), while the notations 
a\ and P\ will still be used where A is the level index. 
Writing the density matrix in terms of the positive-P dis¬ 
tribution (Eq. [2]) transforms the master equation Eq. (pQ) 
into a Fokker-Planck equation: 


dV{a) = f d(a A V(a)) dijh (P(a)) \ 

dt 2r “ V da\ dB\ ) 
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( 4 ) 


According to the Ito calculus, a Fokker-Planck equation 
of the form 

dV(a) _ _ y-^ d ( f n (d)V(d )) + y-^ d 2 (M nrn (a)V(a)) 
dt ^ da n ^ da n ddm 

n nm 

(5) 

with the symmetric matrix M(a) = ^B(o)B(df) T is 
equivalent to the set of Langevin equations [27j 


= fn(di) + B(d) nrn r] rn (6) 

where rj m are independent stochastic Gaussian noise 
terms, defined by (r]m(t)Vn(t')) = — t'). Observ¬ 

able quantities are obtained from averaging the Langevin 
equation over a stochastically generated ensemble. In 
particular, the average occupations (d A a A ) = (n\) and 
second order correlations (d A d A a A a A ) are given by: 


(«a«a) = > ( 7a ) 

(a{a[a x a x ) = {P\P\a\a\) ■ (7b) 


The normalized second order correlation function is then 
defined as g 2 , aa = (at^a A a A )/(n A ) 2 . % solving Eq. (gj) 
numerically, we are able to study the quantum statistical 
properties of bosonic cascades under conditions of macro¬ 
scopic occupations, with several millions of particles. 

Two-Level Case .— Considering the simplest case of 
two levels, we obtain the time dependence of the average 
populations and second order correlation functions shown 
in Figs. [2^ and b, respectively, for the case where the 
upper state is initially a coherent state populated with 
( 77 - 2 ) = 5 x 10 7 particles and the ground state is in the vac¬ 
uum. The dynamics of relaxation can be well understood: 
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FIG. 2. Evolution of two bosonic levels following a co¬ 
herent pulsed excitation. a) Average level occupations, 
(n\(t)). The dashed curve shows the exponential decay of 
the total population, b) Second-order correlation function, 
The dashed curve shows the cross-correlation function 
(a^a^aifa) / ((ni) (n 2 )) • A coherent state of ( 712 ) = 5 x 10 7 
particles is assumed in the highest level and we set Vr — 
8.3 x lO -7 H. 


the population from the upper level decreases in time due 
to radiative emission and transfers to the lower level. At 
early times, the upper state populates the ground state 
in a thermal (or chaotic) quantum state, since it merely 
provides an incoherent input that increases the average 
population without developing any other independent ob¬ 
servable [20]. Consequently, the particle statistics of the 
lower level is initially g2,\=i(t) = 2 which is the well 
known second-order correlation for an incoherent gas of 
bosons, due to their indistinguishability. Then, as popu¬ 
lation increases, stimulated emission sets in and the up¬ 
per state now empties much faster and predominantly 
into the ground state. This results in the rapid growth 
of the ground state population until the upper state is so 
much depleted that it cannot compensate for the ground 
state’s radiative losses, at which stage the ground state 
starts to decay through its own radiative channel. In 
this buildup phase of the ground state, coherence also 
grows, as seen in the transition from #2,a=i(£) = 2 to 
# 2 ,a=i(£) = 1- Note that the state remains coherent from 
there on as radiative decay alone does not dephase the 
condensate. More striking is the statistics of the upper 
state. While its second-order correlation function was 
initially unity, as befits a coherent state, and remained 
essentially unaffected in the first phase of radiative decay 
and spontaneous emission into the ground state, there is 
a pronounced super-bunching g 2 ,\= 2 (t) 2 that forms 

when the state is rapidly emptied. Admirably, just as the 
statistics of the ground state remains equal to one until 
complete evaporation, the upper state’s statistics also re¬ 
mains pinned at the value it reached when it completed 
its transfer. This well defined and high value of $2 for a 
state asymptotically approaching the vacuum is due to 
the mathematical limit of the two vanishing quantities 
in Eqs. © that exists even for arbitrarily small proba¬ 
bilities of occupation. In practice, numerical simulations 


are less stable in this region and experimental measure¬ 
ments would also be increasingly difficult. Nevertheless, 
the most interesting phenomenology is the superbunching 
that develops in the phase where the upper state quickly 
empties into the ground state. In this process, where one 
condensate is sucked by another one, the photons emit¬ 
ted radiatively by the upper state will indeed exhibit a 
physically observable super-bunched statistics. 

The reason for this peculiar behaviour is to be found 
in the mechanism of coherence buildup. Equation (pQ) 
is equivalent to a quantum Boltzmann master equa¬ 
tion 18, 1 28], in which formalism coherence—as measured 
by Glauber’s correlation functions—builds up thanks to 
population correlations developed by the dynamics, even 
in the absence of interactions or external potentials. In 
the conventional case of Bose condensation, the Pois¬ 
son fluctuations of a single condensed state (usually the 
ground state), which leads to g n = 1 for all n, are due 
to this single state adquiring the fluctuations of a macro¬ 
scopic system [2j|. Since the central limit theorem states 
that the distribution of a large number of random vari¬ 
ables (the excited states) is a Poisson distribution, so 
does also fluctuate the condensate in one single state (the 
ground state). Each excited state taken in isolation con¬ 
tributes very little to the condensate in a macroscopic 
system and its statistical properties are not significantly 
altered. In bosonic cascades, however, the asymmetry 
between ground and excited states is lost since there can 
be a few macroscopically occupied excited states. This 
is a peculiar configuration where coherence is adquired 
from another single coherent state, rather than from a 
macroscopic ensemble of weakly occupied thermal states 
(the statistical properties of which do not matter, still 
following the central limit theorem, as long as they obey 
generic conditions of independence and normalization). 
This peculiarity is the reason why the excited state de¬ 
velops such a pronounced superbunching when acting as 
a reservoir for another condensate to grow in another sin¬ 
gle state. The fast loss of coherence from a single state to 
provide for the coherence of another single state, results 
in the superbunching, or extreme chaos, for the provider 
that substitutes a macroscopic environment. 

Further consequences of this mechanism are that the 
superbunching should become stronger for bigger initial 
populations and weaker for more states in the cascade. 
This is expected from the greater dissimilitude of the 
bosonic cascade from a single condensate and a macro¬ 
scopic reservoir: the more levels there are and/or the 
less they are occupied, the more the system resembles 
the usual scenario. Indeed, numerical simulations show 
that this is the case. The dependence of the second or¬ 
der correlation functions, together with the occupation 
numbers, on the initial occupation number is illustrated 
in Fig. [3] and confirms the expected behaviour of # 2,22 
with higher populations. The lowest mode maintains 
the behaviour of being initially thermal, but then be- 
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FIG. 3. Power dependence of the two-level cascade, a) Av¬ 
erage occupation of the lowest level, b) Average occupation 
of the highest level, c) < 72,11 of the lowest level, d) # 2,22 of 
the highest level. Contours correspond to the values shown 
in the colour bars. The horizontal dotted line corresponds to 
the initial occupation considered in Fig. [2] In (d) the value 
of # 2,22 is not defined when ( 712 ) is small. 


comes a coherent state upon being highly populated, as 
the upper state makes the transition from coherent to 
super-bunched light. The greater and sooner the super- 
bunching, the higher the initial population. The effect of 
adding more levels in the cascades also results in weaker 
super-bunching, as expected, but as it also comes with 
notable features of its own, we will discuss it separately. 

Multi-Level Case. — We now consider a larger cascade 
composed of five levels, still with a macroscopic occu¬ 
pation of the highest excited state ((n$(t = 0)) = 10 8 ) 
and all others empty. The results are shown in Fig. [H 
Looking first at the populations (Fig. 0ki), we see that 
there is a staggered oscillation in the successive levels, 
corresponding to the transfers of particles down through 
each level of the cascade. The beginning of the relax¬ 
ation, between the two-highest excited states, is similar 
to the case already studied, but since the recipient for 
the highest excited state is also the source for another 
state below, the dynamics is echoed down the ladder un¬ 
til it reaches the final ground state. In the process, the 
dynamics slows down and is tamed in intensity, as parti¬ 
cles are constantly lost to the environment. The second 
order correlation functions also exhibit a staggering ef¬ 
fect, which can be understood as a generalization of the 
behaviour observed in Fig. [2] for the two-level case, but 
with a richer phenomenology. The super-bunching is still 
observed when the feeding condensate is sucked into the 



FIG. 4. Evolution of five bosonic levels following a coher¬ 
ent pulsed excitation, a) Average level occupations, ( n\{t )). 
The dashed curve shows the exponential decay of the total 
population, b) Second-order correlation function, g 2 ,x(t). 


growing one but this process gets disturbed when the 
latter condensate becomes in its turn the feeding singly- 
occupied macroscopic state that trades its coherence to 
the state below it. This releases the drain on the former 
condensate, which therefore relaxes the super-bunching 
of what is left of the particles there. In the phase of con¬ 
densation, the state that grows its coherence does so not 
only at the expense of its exciting state, but also of the 
state below. This results in a smaller plateau of super- 
bunching for the latter state that sandwiches the conden¬ 
sate together with the plateau of large super-bunching 
from the exciting state. As a result, the particle statistics 
for each state is an intricate sequence of several plateaus 
joined by abrupt jumps as the condensate starts to form 
or starts to empty. For a generic state—that is, one that 
is not too close from the most excited state or from the 
ground state—the statistical relief has five plateaus: i) a 
thermal state starting with, and growing from, the vac¬ 
uum before the cascade is started; ii) a plateau of “small 
super-bunching” as it provides coherence from below to 
its exciting state; iii) a coherent plateau as the state is 
building its own coherence from the excited state now 
in its process of avalanche; iv) a plateau of “big super- 
bunching” as the state feeds the state below, becoming 
super-chaotic in the transfer of coherence; finally, v) a 
thermal state as the process got transported to states 
below, with the next state now in its phase of “big super- 
bunching” and that two states below in that of building 
its coherence. The most relevant plateau is that of stage 
iv), as it is present for all the excited states and has the 
stronger signal. All the excited states except that im¬ 
mediately above the ground state ultimately revert to 
a thermal state, once the peculiar dynamics of bosonic 
cascade is long gone. The complex concordance of these 
several stages in statistics with the corresponding ones 
in populations can be observed in Fig. [ 4 ) 3 , showing the 
beautiful and peculiar dynamics of quantum statistics in 
bosonic cascades. 
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Conclusion .— By solving the Fokker-Planck equation 
for the positive-P representation of a ladder of bosonic 
states coupled as nearest-neighbors, we could compute 
the quantum statistics of a bosonic cascade with sev¬ 
eral levels and for macroscopic occupations. We have 
shown how the peculiar nature of this system, where co¬ 
herence is exchanged between single states, results in a 
super-bunching in the phase where one condensate emp¬ 
ties suddenly into the state below. The observation of 
this strong and striking quantum optical effect could be 
made in wide parabolic quantum wells where excitons 
are confined as whole particles. It will not only illus¬ 
trate fundamental features of coherence buildup in non¬ 
interacting bosonic gases, but also provide new venues for 
imaging applications, in addition to the prospects of im¬ 
plementations in the terahertz bandwidth already offered 
by these systems. 
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